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摘 要 : 空气 动力 学 粗糙 度 是 衡量 地 球 表 


与 大 气 之 间 动 量 和 能 量 交换 的 重要 参数 ,对 于 研究 各 


种 地 表 过 程 和 气候 变化 至 关 重 要 。 琐 感 技术 作为 远 距 离 监测 手段 ,研究 空气 动力 学 粗糙 度 时 其 优 
势 在 于 高 时 效 、 高 经 济 效益 ,能 实现 区 域 或 大 空间 尺度 的 动态 监测 ,因此 利用 还 感 技术 估算 空气 动 
力学 粗糙 度 成 为 热点 问题 。 通 过 系统 阐述 近年 来 国内 外 空气 动力 粗糙 度 研 究 进展 ,重点 介绍 了 利 
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空气 动力 学 粗 糖度 的 方法 ,对 各 种 估算 方法 的 优势 和 不 足 进行 了 总 


结 , 分 析 了 气象 因素 和 地 表 粗 糙 元 形态 特征 因素 对 空气 动力 学 粗糙 度 的 影响 ,进而 对 还 感 技术 在 
该 领域 的 应 用 做 出 展望 , 旨 在 为 空气 动力 学 粗糙 度 遥 感 监测 的 研究 提供 思路 。 
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空气 动力 学 粗糙 度 指 近 地 表 大 气 在 中 性 稳定 
的 条 件 下 ,风速 廓 线 上 风速 为 零 的 几何 高 度 。 由 于 
地 表 粗 燃 元 素 的 存在 ,计算 对 数 风速 剖面 的 有 效 参 


下 垫 面 的 空气 动力 学 粗糙 度 设 为 定 值 ,忽略 了 特定 
土地 利用 类 型 固有 的 时 空 变化 “”。 通 感 数据 覆盖 
范围 大 、 观 测 周 期 短 且 有 多 种 途径 获得 ,为 研究 大 


考 高 度 向 上 移动 至 某 一 高 程 ,这 个 高 度 被 定义 为 零 
平面 位 移 高 度 , 它 反映 了 粗糙 元 素 对 地 表 动 量 吸收 
的 特征 。 空 气动 力学 粗糙 度 表征 地 表 不 同 下 垫 面 
与 近 地 表 气流 的 相互 作用 和 物质 、 能 量 交换 ,对 
土壤 风蚀 有 直接 影响 ,是 研究 干旱 和 半 干 旱地 区 水 
土 流失 的 关键 参数 ” ,也 广泛 应 用 于 区 域 地 表 通 量 
反 演 .化 学 输 运 模型 和 生态 系统 模型 中 “"。 因 此 ， 
空气 动力 学 粗糙 度 的 参数 化 是 值得 深入 人 研究 的 
课题 。 

目前 , 随 着 国内 外 学 者 对 地 表 空 气动 力学 领域 
的 深入 开展 ,估算 区 域 尺 度 上 非 均 匀 地 表 的 空气 动 
力学 粗糙 度 成 为 热点 问题 。 基 于 实测 数据 的 估算 
方法 是 验证 遥感 方法 准确 性 的 前 提 *" ,主要 有 和 冠 
层 高 度 固 定 比例 法 、 风 洞 试验 法 和 野外 实验 法 。 后 
来 发 展 的 查 表 法 能 实现 更 大 范围 空气 动力 学 粗粮 
度 的 估算 ,该 方法 基于 土地 利用 类 型 ,将 同一 性 质 
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空间 尺度 的 空气 动力 学 粗糙 度 的 动态 变化 开辟 道 
路 。 近 20 a 里 ,国内 外 学 者 对 遥感 技术 应 用 于 空气 
动力 学 粗糙 度 的 参数 化 进行 了 多 方面 研究 ””” 。 
本 文 概述 了 利用 遥感 技术 估算 植被 下 垫 面 空气 动 
力学 粗 燃 度 的 方法 。 


1 基于 实测 数据 的 方法 


1.1 基于 冠 层 高 度 固定 比例 法 

用 植被 平均 高 度 佑 算 空 气动 力学 粗糙 度 是 最 
常用 的 方法 之 一 ,对 于 植被 下 热 面 ,植被 高 度 是 地 
表 粗 燃 状 况 及 气流 拖 趾 力 大 小 的 影响 因子 。 
在 相关 研究 中 ,按照 简单 的 经 验 法 则 将 植被 高 度 与 
空气 动力 学 粗糙 度 联 系 起 来 近似 计算 hh, 
Yang 等 ”对 具有 复杂 几何 和 空间 结构 的 北方 森林 
冠 层 的 空气 动力 学 粗糙 度 进行 参数 化 ,使 用 植被 高 
度 ( 刀 ) 来 简化 表示 空气 动力 学 粗糙 度 (a) ,计算 公式 


基金 项 目 : 国家 自然 科学 基金 面 上 项 目 (41871003 ) ;烟台 市 科技 创新 发 展 计划 项 目 (2022MSGY062 ) 资 助 
作者 简介 : TE FEE (1998) , 女 ,硕士 研究 生 ,主要 从 事 草原 植被 遥感 监测 研究 .E-mail: lxywork1@163.com 
通讯 作者 : 王 静 珊 (1987- ) , 女 ,博士 ,副教授 , 主要 从 事 植 被 遥感 监测 .沙漠 化 过 程 研究 . E-mail: wangjp@ldu.edu.cn 


202304.00892v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


408 T 4 wA 4635 


如 下 : 
z0=0.125H (1) 

Masseroni 等 汪 对 稻田 在 整个 农业 季节 的 零 平 
面 位 移 高 度 和 空气 动力 学 粗糙 度 进 行 测量 发 现 , 在 
整个 生长 季 ,2 个 参数 与 三 有 良好 的 相关 性 (R?> 
0.9) ; Santos 等 站 基 于 卫星 的 能 量 平衡 模型 (MET- 
RIC ) 研 究 橄榄 园 的 蒸 散 量 ,模型 中 空气 动力 学 粗糙 
度 (z) 的 固定 函数 为 : 

Zza =0.068H (2) 

Van der Graaf 等 “在 用 于 研究 欧洲 西北 部 活性 
扼 的 地 表 - 大 气 交 换 模型 中 ,空气 动力 学 粗糙 度 作 
为 模型 输入 参数 ,其 大 小 取决 于 植被 高 度 。 
12 野外 实验 法 

由 于 观测 条 件 和 实验 观测 仪器 的 差异 ,在 野外 
实测 气象 数据 时 主要 发 展 了 风速 廊 线 法 、 涡 动 协 方 
差 (EC ) 法 和 大 孔径 闪烁 仪 (LAS) 法 。 利 用 野外 实测 
数据 的 推算 方法 有 最 小 二 乘 拟 合 迭代 法 AF BG 
法 温度 方差 法 (TVM) 和 Martano 法 1 
1.2.1 Mik BRE 气象 站 中 常用 到 的 监测 设备 主 
BOA UTR AG ait AUT FEI A . 温 湿 度 传感器 等 , 传 
感 器 将 对 应 气象 要 素 值 存储 在 数据 采集 器 中 。 利 
用 气象 观测 资料 计算 空气 动力 学 粗糙 度 时 ,采用 
Monin-Obukhov 相似 度 理论 ,基于 不 同 高 度 同一 时 
刻 测 量 的 风速 剖面 和 温度 剖面 数据 ,根据 微 气象 学 
原理 及 其 一 系列 计算 公式 ,通过 最 小 二 乘法 拟 合 迭 
TRAE BIAS BH J) 2f Ub BE Eee DUREE nmm. Hm, 
对 风速 廊 线 的 气象 剖面 时 间 序 列 数 据 进 行 处 理 , 在 
中 性 稳定 大 气 条 件 下 ,计算 公式 如 下 : 


| (3) 
SOP iu UR (m s );w: 为 摩 阻 速度 (ms ) 5k H 
Von Karman 常数 ( 取 值 0.4);a 为 空气 动力 学 粗糙 度 
(m) ;d 为 零 平 面 位 移 m);z 为 参考 高 度 m)。 对 测 
定 结果 用 最 小 二 乘法 进行 回归 分 析 ,确定 最 佳 拟 合 
关系 式 来 求 出 z。 

薛 唱 等 采用 风速 廓 线 仪 采集 不 同 高 度 草原 
灌木 带 的 风速 廓 线 数据 ,采用 最 小 二 乘法 对 数据 进 
行 计算 得 到 空气 动力 学 粗糙 度 ;Yu 等 “利用 干旱 和 
半 湿 润 气候 条 件 下 自动 气象 站 的 风速 资料 ,采用 和 迭 
代 法 计算 2 个 地 区 时 间 序 列 上 的 空气 动力 学 粗糙 度 。 
1.2.2 1% ha Z (EC) & fo K 7042 A HRA (LAS) 
法 EC 和 LAS 法 采用 Monin-Obukhov 相似 度 理论 ， 


从 地 表 通 量 的 角度 出 发 ,得 到 计算 空气 动力 学 粗糙 
度 的 变量 ,结合 当地 气象 数据 ,进而 求 得 空气 动力 
学 粗 烽 度 的 方法 ”™'"。2 种 方法 区 别 在 于 :一 是 
LAS 与 EC 相 比 ,由 于 光路 从 数 百 米 到 数 千 米 不 等 ， 
能 够 测量 大 面积 的 庙 流 通 量 ,更 适合 计算 大 空间 斥 
度 上 的 非 均 匀 下 垫 面 空气 动力 学 粗糙 度 ”; 二 是 在 
计算 通 量 时 ,EC 数据 包括 风速 .大 气 稳 定性 、 摩 擦 速 
度 等 计算 空气 动力 学 粗糙 度 的 变量 ,LAS 数 据 包括 
大 气 稳定 性 和 摩擦 速度 这 2 个 变量 ,而 风速 数据 要 
根据 气象 站 观测 来 获得 。 

EC 观测 系统 是 一 种 微 气象 学 的 测量 方法 ,采用 
EC 原理 ,利用 快速 响应 的 传感器 来 测量 大 气 与 下 垫 
面 间 的 物质 交换 和 能 量 交 换 , 它 是 目前 地 表 通 量 观 
测 精 度 最 高 且 最 常用 的 仪器 。EC 通 量 观测 站 可 以 
观测 到 瞬时 风向 、 瞬 时 风速 ,温度 、 水 汽 感 热 通 量 和 
痕 量 气体 浓度 数据 ,是 一 种 直接 测算 通 量 的 标准 方 
法 。LAS 由 发 射 仪 和 接收 仪 2 部 分 组 成 ,接收 仪 接 
收 到 光 程 上 大 气 波动 影响 的 发 射 波束 ,并 用 折射 指 
数 的 结构 参数 来 表达 大 气 的 消 流 强度 ,进而 推算 滑 
流通 量 。 当 光束 在 扰动 大 气 中 沿 直线 传播 时 ,受到 
扰动 引起 传播 光束 强度 的 波动 ,通过 接收 到 光 程 上 
大 气 波动 光束 ,并 用 折射 指数 结构 参数 ( C7 ) 来 反映 
大 气 的 济 流 强度 ”1。 

根据 Monin-Obukhov 相似 度 理论 ,Sun 等 ”利用 
EC 数据 和 LAS 数 据 计算 了 青藏 高 原 北部 的 空气 动 
力学 粗 焰 度 , 空 气动 力学 粗 燃 度 的 计算 公式 为 : 


eu ed 
2 i «v L ) (4) 


式 中 :za 为 空气 动力 学 粗糙 度 Am) 5 ul) H EC EÈ LAS 
测量 高 度 z 处 的 水 平 风速 (ms 0 5d 为 零 平 面 位 移 
(m) ,设置 为 植被 高 度 的 2/3;L 为 Monin-Obokhov 长 
RE (m) ;ww 为 摩 阻 速度 (ms ') ;名 ,为 大 气 稳定 性 修正 
PRA sk Von Karman 常数 , 取 值 0.4。 

Bg se re Sel A A A eA A aA 
资料 估算 农田 .草地 下 热 面 的 空气 动力 学 粗糙 度 时 
空 变化 特征 ,并 将 估算 值 与 其 他 参数 化 方案 进行 比 
较 ;Liu 等 ”利用 EC 法 测量 值 计算 的 空气 动力 学 粗 
Jis HE CR: ur Zs Uy pa 989 A A ERN. RE 28 ^62] 73 
学 粗糙 度 ;Zhang 等 3 利用 三 维 超声 波 风速 仪 获取 
中 国 北方 风蚀 区 5 个 站 点 的 近 表面 风速 ,计算 草地 、 
农田 流动 沙 地 和 戈壁 4 种 地 表 类 型 的 空气 动力 学 
粗糙 度 ,发 现 草地 的 地 表 粗 糙 度 最 大 ,其 次 是 农田 、 
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沙 地 和 戈壁 。 相似 性 原则 ,有 别 于 地 表 真 实 状况 ,研究 发 现 风 洞 


1.3 风 洞 试验 法 

风 洞 试验 通过 在 实验 室 模 拟 自 然 风 对 不 同形 
状 及 植被 覆盖 密度 情况 下 土壤 的 风蚀 作用 ,得 到 不 
同 高 度 的 风速 ,再 由 中 性 稳定 大 气 条 件 下 的 计算 公 
式 [ 式 (3)] 得 到 空气 动力 学 粗 烟 度 "3 。 风 洞 
分 为 过 渡 段 整流 段 收缩 段 和 试验 段 ,同时 配套 有 
风速 测试 控制 系统 ,其 中 试验 段 模拟 自然 风 , 能 六 
生 自 由 旋涡 气流 和 稳定 流动 的 气流 场 "。 风 洞 试 
验 设 置 模型 作物 的 高 度 abis RE 、 株 行距 等 重要 参 
数 ,模拟 地 表 植 被 不 同 粗糙 程度 下 对 风 的 阻挡 作 
用 ,研究 空气 动力 学 粗糙 度 时 与 冠 层 高 度 固定 比例 
法 相 比 更 为 准确 。 

植被 作为 主要 的 风蚀 因子 之 一 ,近年 来 国内 外 
学 者 开展 了 风 洞 模拟 实验 用 于 研究 植被 表面 粗糙 
度 对 土壤 风蚀 作用 的 影响 。 张 春来 等 利用 风 洞 
试验 研究 了 耕作 土壤 表面 的 空气 动力 学 粗糙 度 ; Y 
日 娜 等 中 基于 风 洞 模拟 试验 ,研究 沙漠 地 区 油 沙 豆 
种 植 示范 区 在 不 同 留 茬 模式 下 近 地 表 风速 变化 及 
防风 效能 ,采用 风速 廓 线 法 计算 空气 动力 学 粗 
Ts BE 

移动 式 风蚀 风 洞 的 出 现 解 决 了 传统 风 洞 试验 
法 取样 时 对 原 地 表土 造成 的 扰动 的 问题 。 刘 艳 萍 
等 3 利用 移动 式 风蚀 风 洞 研究 草原 区 植被 对 土壤 
风蚀 的 影响 ; 赵 永 来 等 利用 移动 式 风 蚀 风 洞 分 析 
了 燕麦 残 茬 下 垫 面 的 空气 动力 学 粗糙 度 的 变化 。 
1.4 方法 评述 

冠 层 高 度 固定 比例 法 采用 单一 变量 的 简化 经 
验 关 系 式 ,与 植被 平均 高 度 相关 ,但 相关 人 研究 表明 ， 
气象 因子 和 地 形 起 伏 等 对 植被 下 热 面 空气 动力 学 
粗糙 度 均 存 在 影响 ?2 , 且 该 方法 一 般 适 用 于 农 
业 地 区 的 均匀 植被 地 表 ,而 区 域 或 大 空间 尺度 的 植 
被 具有 空间 异 质 性 ,因此 在 研究 空气 动力 学 粗糙 度 
时 存在 一 定局 限 性 。 野 外 实验 法 对 数据 质量 要 求 
严格 ,在 野外 设置 不 同 站 点 获取 气象 数据 ,一 般 用 
来 作为 地 面 验证 的 真 值 ,但 自然 条 件 和 仪器 系统 的 
微小 影响 会 使 其 计算 结果 产生 较 大 误差 , 且 结 果 只 
能 代表 一 定 范围 源 区 ,不 能 推广 到 大 空间 尺度 的 应 
用 。 风 洞 试验 法 考虑 了 气流 与 地 表 粗 糙 元 对 空气 
动力 学 粗 烟 度 的 综合 作用 ,模拟 理想 状态 下 不 同 风 
速 和 地 表 和 覆盖 状况 下 的 风蚀 作用 ,进而 研究 空气 动 
力学 粗糙 度 ,但 不 足 之 处 在 于 风 洞 试验 是 采取 流 场 


实验 与 现场 实测 的 风沙 结果 相似 性 较 低 ”, 且 所 需 
仪器 设备 成 本 较 高 且 耗 时 耗 力 ”。 


2 遥感 方法 


遥感 影像 覆盖 范围 大 、 观 测 周期 短 上 且 有 多 种 途 
径 获 得 ,因此 在 研究 大 空间 尺度 的 空气 动力 学 粗糙 
度 时 具备 突出 的 科学 性 和 技术 优点 。 遥 感 方法 通 
过 反 演 地 表 参 数 或 者 大 气 参数 输入 到 模型 计算 空 
气动 力学 粗糙 度 ' ?0 。 

2.1 植被 指数 模型 法 

空气 动力 学 粗糙 度 与 下 垫 面 粗糙 元 的 地 表 特 
征 有 关 , 植 被 与 水 体 、 裸 地 等 下 垫 面 相 比 ,形态 特征 
复杂 .动态 变化 明显 ,因此 在 遥感 中 引入 植被 指数 
估算 空气 动力 学 粗糙 度 是 非常 有 必要 的 。 植 被 指 
数 是 根据 植被 的 光谱 特性 ,基于 多 种 波段 组 合 ,对 
地 表 植 被 状况 的 度量 。 基 于 光学 卫星 数据 的 反 演 
方法 一 般 是 建立 植被 指数 与 空气 动力 学 粗糙 度 E 
平面 位 移 高 度 的 关系 模型 ,再 通过 经 验 模型 转换 为 
空气 动力 学 粗糙 度 。 常 用 的 植被 指数 有 归 一 化 植 
被 指数 (NDVI) 、 叶 面积 指数 (LAT) .锋面 面积 指数 
(FAI). 

2.1.1 归 一 化 植被 指数 (NDVI) NDVI 是 植物 生长 
状态 以 及 植被 空间 分 布 密度 的 最 佳 指示 因子 。ND- 
VI 是 根据 近 红 外 波段 反射 率 (psi) 和 红 光 波段 反射 
率 (pii) 计 算出 来 的 ,计算 公式 如 下 : 
NDVI- Ps Pea (5) 
Pair TP red 

利用 遥感 技术 研究 植被 下 垫 面 的 空气 动力 学 
粗糙 度 时 , NDVI 的 应 用 较为 广泛 >”** EE 
AEU H Landsat TM 资料 计算 NDVI, 将 其 与 黑河 试验 
点 的 潮流 通 量 观测 所 估算 的 空气 动力 学 粗糙 度 建 
立 经 验 关系 式 ;Gupta 等 ”使 用 NOAA-AVHRR 数据 
计算 NDVT, 结 合 地 表 过 程 参数 建立 印度 半 王 旱地 
区 空气 动力 学 粗糙 度 与 NDVI 之 间 的 关系 式 ;Xing 
等 呈 对 比 不 同 植被 指数 与 空气 动力 学 粗糙 度 的 关 
系 ,发现 NDVI 是 草地 空气 动力 学 粗糙 度 的 敏感 指 
Ek; Abbas 等 5 得 到 了 描述 NDVI 与 空气 动力 学 粗糙 
度 关系 的 最 佳 最 终 数 学 模型 ,用 于 绘制 整个 伊拉克 
地 区 的 粗糙 度 长 度 图 中。 

2.1.2 叶 面 积 指 数 (LAT) LAI 量 化 了 陆地 植被 的 
绿叶 面积 ,是 植被 冠 层 结构 和 功能 的 基本 属性 ””。 
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大 量 研究 表明 ,空气 动力 学 粗糙 度 与 植被 冠 层 结构 
参数 有 关 ,因此 在 反 演 空气 动力 学 粗糙 度 时 ,LAI 在 
一 定 程度 上 比 NDVI 具 有 更 强 的 相关 性 , 零 平 面 位 移 
高 度 和 空气 动力 学 粗糙 度 是 LAI 的 强 冰 数 ?935 "5 。 
张杰 等 所 根据 空气 动力 学 粗糙 度 与 风速 .植被 密度 
和 LAI 的 关系 ,采用 遥感 方法 反 滨 LAI 和 风速 ,估算 
稀 玻 植被 区 空气 动力 学 粗糙 度 ; Cho 57 HF] LAI 与 
冠 层 高 度 的 比值 表征 冠 层 结构 的 复杂 性 ,来 检验 空 
气动 力学 粗糙 度 与 反照 率 的 关系 ;Alekseychik 45: ^^ 
将 空气 动力 学 粗糙 度 与 LAI 表 示 为 线性 函数 关系 ; 
Sun 等 3 利用 前 人 开发 的 经 验 模型 和 LAI 数 据 , 研 究 
青藏 高 原 北部 小 空间 尺度 上 空气 动力 学 粗糙 度 的 
空间 分 布 ;Liu 等 利用 满 流 传输 模型 ,根据 LAI 和 
植物 功能 类 型 有 关 的 冠 层 形态 特征 ,计算 了 1982 一 
2017 年 的 全 球 零 平面 位 移 高 度 和 空气 动力 学 粗糙 
度 这 2 个 空气 动力 学 参数 。 
2.1.3 锋面 面积 指数 (FAI) FAI 表 示 为 粗糙 度 元 素 
的 锋面 积 与 单位 地 面 面 积 的 比值 ,后 来 用 冠 层 面积 
BACCA ) 或 单位 地 面 面 积 中 所 有 冠 层 元 素 的 总 单 面 
面积 来 替代 FAI,4 包 括 所 有 吸收 动量 的 冠 层 元 素 。 
Jasinski 等 ”建立 了 自然 植物 或 随机 分 布 植物 的 
FAI 和 4 的 参数 化 方法 ,并 将 其 应 用 于 利用 卫星 网 像 
估计 局 部 空气 动力 学 粗糙 度 ;Schaudt 等 所 采用 Rau- 
pach 提出 的 空气 动力 学 粗糙 度 与 FAI 的 函数 关系 
式 , 利 用 卫星 数据 反 演 空气 动力 学 粗糙 度 , 当 FAI= 
0.152 时 ,表示 为 : 
= = 5.86 exp(-10.94'")a'*+0.00086 — (6) 


M FATDO0.152 时 ,表示 为 : 


Zo = 0.093711 0 —exp(—10.94''")] +0.00368 (7) 


式 中 :za 为 空气 动力 学 粗糙 度 (m) ;h 为 植被 高 度 
(m); 和 为 单位 面积 的 锋面 面积 指数 。 

杨 阿 强 等 “采用 前 人 改进 的 形态 学 模型 ,根据 

4 与 LAI 之 间 的 关系 ,采用 遥感 数据 源 对 中 国 东部 

的 空气 动力 学 粗糙 度 进 行 估算 ,并 人 研究 其 时 空 变 

化 ,其 中 A 定义 为 单位 面积 上 植被 对 风 的 阻挡 面积 ， 

A=LAI+Lg+Ls (8) 

式 中 : Le 为 植株 地 上 茎 面积 指数 ; Ls AHS AAR 

指数 。 

利用 LAI 和 经 验 参 数 计算 Lg 和 Ls 二 者 的 和 ，, L 

为 二 者 的 和 ,其 中 为 第 n 个 月 , o 为 不 同 植被 类 型 的 


固定 参数 ,1- a 为 每 月 枯 叶 凋落 比率 ,计算 公式 如 下 : 
L'-lg*ls- 
(9) 


214 双向 反射 率 分 布 函 数 指数 从 单一 观测 角度 
计算 的 植被 指数 在 反 演 与 空气 动力 学 粗糙 度 密切 
相关 的 三 维 植被 结构 方面 能 力 有 限 ,一 个 潜在 的 解 
决 方案 是 引入 多 角度 光学 遥感 ,标准 MODIS 
BRDF/ 反 照 率 (MCD43) 实 现 了 半 经 验 BRDF 模 型 数 
据 产品 。Xing 等 ”将 草地 下 扑面 的 BRDF_R 指数 与 
空气 动力 学 粗糙 度 的 关系 表示 为 指数 过 函数 的 
XX GR D ,其 中 Ross-Li 模 型 采用 3 个 核 的 线性 加 权 
和 形式 ,将 二 向 反射 率 表示 为 3 项 的 总 和 : 

Pp(0.,0,,9)=k +k F (0,0,0) - k,F,(6,.0,, 9) (10) 
式 中 : p(0.,0.,9) 为 二 回 反 射 率 ; 9. 为 太阳 天 顶 角 ; 0, 
为 观测 天 顶 角 ; o 为 太阳 和 观测 点 的 相对 方位 角 ; 
所 和, 分别 为 体 散 射 核 和 几何 光学 散射 核 , 均 是 
0,.0, p 的 函数 ; ky ki 和 ,为 各 个 核 的 系数 ,分 
别 表示 各 向 同性 散射 、 体 散射 .几何 光学 散射 所 占 
的 权重 。 基 于 式 (10) ,反映 地 表 几 何 粗 糙 度 的 近 红 
外 双向 反射 率 分 布 玫 数 指数 (BRDF_R) ,计算 公式 
如 下 : 


NS 


BRDF_R=k,/k, (11) 

通过 遥感 数据 提取 植被 指数 进行 空气 动力 学 
粗糙 度 的 研究 不 断 更 新 ,模型 形式 多 样 ,近年 来 研 
究 方 法 归纳 总 结 见 表 1。 
2.2 激光 雷达 方法 

空气 动力 学 粗糙 度 受 地 表 粗 糙 度 元 素 的 高 度 、 
几何 形状 、 密 度 等 影响 ,这些 元 素 包 括 植被 以 及 微 
WAZ NUL HUE RAE 2005055 激光 雷 达 在 描 
述 地 表 和 粗糙 元 特征 方面 显示 出 巨大 潜力 ,在 
研究 空气 动力 学 粗 烟 度 时 比 光 学 卫星 具有 一 定 技 
术 优 势 , 按 搭载 平台 来 看 主要 有 机 载 激 光 雷 达 和 卫 
星 雷达 。 机 载 激 光 雷 达 是 一 种 主动 测量 方式 ,利用 
激光 进行 探测 和 测量 ,通过 点 云 数据 获得 陆地 表面 
和 植被 的 三 维 结构 ,为 估算 空气 动力 学 粗糙 度 提 供 
了 新 方法 。Greeley 等 利用 校准 后 的 雷达 后 向 散 
射 系数 与 空气 动力 学 粗糙 度 建立 关系 ,发 现 4 种 类 
型 下 垫 面 均 具 有 很 好 的 相关 性 ,证 明 可 以 利用 雷达 
数据 评估 空气 动力 学 粗糙 度 。 

Colin 45 ? 通过 激光 雷达 测量 的 点 云 数 据 获得 
黑河 流域 的 数字 高 程 模型 (DEM) 和 数字 地 表 模 型 
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表 1 常用 的 空气 动力 学 粗糙 度 参 数 化 方法 


Tab.1 Common parameterization methods of aerodynamic roughness 


指数 植被 类 型 表达 式 参考 文献 
NDVI 农田 zo=exp(-5.5+5.8NDVD) Gupta 等 中 
林地 ,草地 2,70.0206e v Abbas 等 中 
草地 2,-0.0203NDVT7 Xing 等 中 
玉米 .小麦 zo(x, y)=exp[7.13+9.33NDVI(x, y)] 贾 立 等 " 
春玉 米 zo=0.2255NDVI+0.0087 Yu Al 
冬小麦 z0=0.2476NDVI+0.0615 
夏 玉 米 z0=0.2858NDVI+0.1017 
RAGE fet ze MEM Van der Graaf 等 四 
LA 林地 270.3299 1p 52.1713 Schaudt 等 外 
Pp LAL-6.512,-0.17 Alekseychik 4^9 
玉米 zo=z0 +0.3h(CaxLAD™® Lu 等 四 
短 高 山 草 zj-exp[-2.225-0.938/(LAI-0.205)] Sun 等 站 
zo=exp[2.613-0.173/(LAI-0.287)] 
BRDF_R 草地 20=0.0013BRDF_R-0.65 Xing p?" 
FAI 宽 叶 林 、 针 叶 林 z/h-0.0537XF AT" 1.0—-exp(-10.9FAT"**)]-0.00368 Schaudt 4&9 
HDVI 春玉 米 z=0.2236HDVI-0.0279 MEL 
冬小麦 zo=0.2695HDVI+0.0688 
E z=0.2113HDV1+0.0391 
NDVI,u Tit EB 2,-0.53074"—0.3952 Xing 4" 
CI 草地 2,70.0078CI-0.493 Xing $E 
NDVI.a 棉花 zy=exp[0.26(NDVI/a)-2.21] Allen 等 中 


Tu 


注 :za 为 空气 动力 学 粗糙 度 ;NDVI 为 归 一 化 植被 指数 ;Z 为 植物 平均 叶 面 积 指数 ;LAI 为 叶 面积 指数 ;FAI 为 单位 面积 锋面 面积 指数 ;LAI OI P 
草 单一 种 类 的 叶 面 积 指数 ; C, 为 平均 阻力 系数 ;z 为 土壤 表面 粗糙 度 长 
速 


度 沁 为 植被 高 度 ;BRDF_R 为 近 红 外 双向 反射 率 分 布 函 数 指数 刻 为 风 


;HDVI 为 归 一 化 差异 热点 暗 点 指数 ;CI 为 BRDF_R 和 NDVI 结 合 的 综合 指数 ;a 为 反照 率 。 


(DSM) ,将 DEM、DSM 应 用 到 几何 模型 (Macdonald 
模型 .R92 形态 学 模型 ) 和 流体 动力 学 模型 (CFD)， 
发 现 CFD 模型 对 不 同 下 垫 面 的 适用 性 更 强 , 可 以 更 
精细 地 观察 气流 与 粗糙 元 素 之 间 的 相互 作用 ;Li 等 “ 
利用 激光 雷达 和 成 像 光 谱 技 术 获 取 植 物 结构 信息 ， 
估计 以 灌木 为 主 的 早 地 研究 区 域 的 空气 动力 学 粗 
糙 度 ,对 比 了 2 种 参数 化 方法 估计 空气 动力 学 粗糙 
度 的 敏感 性 ; Trepekli 等 :5 利用 无 人 机 载 激 光 雷 达 
获取 植被 表面 的 三 维 结构 数据 ,评估 了 3 种 形态 学 
模型 估算 丹麦 农田 和 周围 常 绿 树木 空气 动力 学 粗 
糖度 的 准确 性 ,其 中 Raupach 模 型 的 估算 结果 为 最 
佳 ,该 方法 基于 植被 几何 形状 和 风 条 件 ,表达 式 
如 下 : 


zo =h(1 - d/h) x exp(-kulu, + yr,) (12) 
dih=1+ {(eof-(2eaar)") e (Ge. ran") (13) 


u./u= min|(c, + c,FAl)”, max(u./u)| (14) 


式 中 :za 为 空气 动力 学 粗糙 度 (m) ; h AE SA E 
TRAE); d 为 零 平 面 位 移 (m) ;为 Von Karman 
常数 ( 取 值 0.4) ; u W UE (m +s") ;ww 为 摩 阻 速度 
(m-s'); vr, ALBEE ESSE WE PR; cj 是 自由 参 
数 ( 取 值 7.5);FAI 由 激光 雷达 获取 的 横 截 面 线 上 的 
正高 度 变化 (Ay ) 除 以 该 截面 长 度 的 距离 (Ax ) 计 
算 ; c, 为 单个 粗糙 度 元 素 的 阻力 系数 ; e, 为 基底 表 
面 在 有 处 的 阻力 系数 。Paul-Limoges 等 5 利用 机 载 
激光 雷达 数据 来 佑 算 冷 杉林 采伐 前 后 的 空气 动力 
学 粗糙 度 ,表明 了 在 复杂 情况 下 使 用 激光 雷达 数据 
估计 空气 动力 学 粗糙 度 的 潜力 。 
2.3 多 源 数 据 结合 方法 

不 同 研究 区 影响 空气 动力 学 粗糙 度 的 主要 因 
RABE ,不 同 遥 感 数据 源 各 具 其 特征 与 优势 , 结 
合 遥 感 多 源 数据 建立 空气 动力 学 模型 显现 其 优势 
所 在 2255 Alen 等 "提出 使 用 反照 率 (w) 有 助 
于 区 分 具有 相似 的 NDVI 但 却 高 度 不 同 的 植被 类 
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型 ,将 NDVI 与 cx 相 结合 ,并 建立 与 空气 动力 学 粗糙 
度 的 关系 式 ( 表 1); 吴 炳 方 等 "综合 考虑 了 植被 信 
息 、 地 形 因子 . 非 植被 覆盖 表面 的 几何 粗糙 度 因 素 
来 估算 空气 动力 学 粗糙 度 ,其 中 植被 指数 NDVI 的 


NDVI 
NDVI rax ~ Aa 2: 


式 中 : zo, 为 植被 区 的 粗 烟 度 (m);NDVI 为 区 域 中 
归 一 化 植被 指数 的 最 大 值 ;NDVL 为 最 小 值 。 

在 地 形 起 伏 的 情况 下 ,考虑 DEM 的 坡度 信息 
(slope) 对 粗糙 度 的 影响 : 


4L =a 1+ =") (16) 


WP: z 为 坡度 因子 的 粗糙 度 (m) ; slope 为 从 数字 

地 形 数据 中 获取 的 坡度 文件 ;a 与 6 为 经 验 系数 。 

从 雷达 影像 中 反 演 得 到 的 表面 几何 粗糙 度 表 
示 非 植被 覆盖 区 域 的 地 表 粗 糙 度 : 

Ig(z,,,) =—1.221 +0.0906(0,) (17) 

式 中 : n 为 非 植被 覆盖 区 域 的 地 表 粗 糙 度 (m) ;or 

为 从 Radarsat 卫星 影像 中 获取 的 后 向 散射 系数 。 区 

域 综合 有 效 粗 糙 度 是 以 上 植被 .地形 和 非 植被 覆盖 

表面 的 几何 粗 烟 度 3 种 因素 的 加 权 函 数 。 

Xing 等 中 发 现 BRDF_R 和 NDVI 都 是 草地 空气 
动力 学 粗 烟 度 的 敏感 指标 ,创新 性 地 将 二 者 结合 开 
发 了 组 合 指数 (CT) ,是 与 生长 季 草 地 空气 动力 学 粗 
烟 度 呈 罕 函数 关系 ( 表 1) ,其 中 ， 

z sqrt(BRDF_R) 
NDVI'° 
Xing 等 又 提出 了 NDVI 和 实测 风速 相 结合 计 

算 空 气动 力学 粗糙 度 的 方案 ( 表 1) ,以 评估 黑河 流 

域 干 旱地 区 绿洲 的 空气 动力 学 粗糙 度 ,该 方案 表现 

出 良好 的 适用 性 ,其 中 ， 

WN=ws™ (19) 

式 中 :WN 为 风速 与 NDVI 的 综合 指数 ;WS 为 NDVI 

数据 日 期 5m 人 处 的 平均 风速 (ms )。 

2.4 方法 评述 
利用 遥感 技术 估算 空气 动力 学 粗糙 度 ,在 区 域 

或 大 空间 尺度 上 具有 明显 优势 ,成 本 低 且 不 会 对 环 

境 产 生 破 坏 。 植 被 指数 模型 将 空气 动力 学 粗糙 度 

与 植被 特性 相 联系 ,但 对 于 非 均 匀 植 被 下 垫 面 ,不 

同和 遥感 产品 得 到 的 同一 植被 指数 可 能 会 有 偏差 H. 

植被 指数 在 非 生 长 季 并 不 适用 ”。 激 光 雷 达 具 有 


z$ 70.014 [os 


CI (18) 


方向 性 好 、 分 辨 率 高 且 抗 干扰 能 力 强 的 技术 优势 ， 
在 获取 地 表 立 体 信息 方面 有 很 大 潜力 ,但 需要 注意 
的 是 ,机 载 激 光 雷 达 会 受 大 气 状况 影响 , 且 获 取 地 表 
植被 结构 信息 会 存在 一 定 的 局 限 性 ,限制 包括 低估 
植被 高 度 和 难以 获得 个 别 粗糙 度 元 素 冠 层 信息 。 
多 源 遥 感 协同 方法 综合 考虑 空气 动力 学 粗糙 度 的 
影响 因子 ,通过 异 质 观测 数据 相 结 合 克 服 单一 数据 
源 存 在 的 缺陷 ,实现 互补 互利 ,来 获得 区 域 最 佳 估 
算 模型 ,但 目前 也 存在 多 源 数 据 在 各 因子 耦合 和 时 
空 分 状 率 匹配 方面 的 问题 。 


3 ”空气 动力 学 粗糙 度 影响 因子 


空气 动力 学 粗糙 度 的 大 小 主要 取决 于 地 表 粗 
糙 度 要 素 的 几何 特征 和 分 布 密度 ,同时 也 受 大 
气 因 子 ( 风 速 `. 风 向、 大 气 稳定 度 等 ) 影 响 。 在 风蚀 
模型 中 ,植被 的 物理 结构 通过 对 气流 施加 阻力 与 风 
相互 作用 ,因此 空气 动力 学 粗糙 度 是 各 种 因子 综合 
作用 的 动力 学 参数 。 
3.1 气象 因素 
311 风速 一 定 地 表 条 件 下 空气 动力 学 粗糙 度 
并 不 是 一 个 常数 ,其 值 随 着 风速 变化 而 有 明显 变 
化 “S331 对 于 植被 覆盖 下 垫 面 , 随 着 风速 的 增 
加 , 风 在 一 定 程度 上 改变 地 表 植 被 的 几何 形状 ,在 
风 的 作用 下 向 同一 方向 倾斜 ,空气 动力 学 粗糙 度 随 
风速 的 增加 而 减 小 达到 一 定 临 界 点 后 稳定 ”52 1。 

张 雅 静 等 ”基于 实测 数据 和 理论 模拟 计算 ,发 
现 不 同 植被 密度 条 件 下 ,空气 动力 学 粗糙 度 均 与 风 
速 时 指数 关系 ;Yu 等 ”在 研究 黑河 流域 不 同时 间 尺 
度 下 空气 动力 学 粗糙 度 的 影响 因素 时 发 现 ,风速 是 
草原 和 森林 山区 下 垫 面 空气 动力 学 粗糙 度 的 主要 
影响 因子 之 一 ,而 对 于 平坦 农田 是 次 要 因素 ;Zhang 
等 3 比较 中 国 北方 风蚀 区 4 种 地 表 类 型 空气 动力 学 
粗糙 度 ,发 现 空气 动力 学 粗糙 度 值 在 低 风速 时 变化 
最 大 ,在 较 高 风速 时 变化 较 小 。 
3.1.2 摩 阻 风速 摩 阴 风速 是 对 由 于 地 面 摩擦 阻力 
存在 而 产生 风速 梯度 的 衡量 ,也 是 气流 对 床 面前 切 
力 的 反映 ,风力 越 强 ,地面 阻力 越 大 。 研 究 发 现 , 摩 
阻 风速 与 平均 风速 存在 线性 关系 “5 。 赵 永 来 等 中 
发 现 各 实验 测 点 的 空气 动力 学 粗糙 度 都 随 着 摩 阻 
风速 的 增 大 而 增 大 ,而 当 风 速 增 大 到 地 表 植 被 倒伏 
时 ,边界 气流 层 切 应 力 增 大 从 而 导致 空气 动力 学 粗 
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糙 度 和 摩 阻 风速 增幅 降低 甚至 变 为 负 值 。 
3.1.3 风向 风向 变化 对 空气 动力 学 粗糙 度 的 影响 
取决 于 地 形 的 不 均匀 性 ,对 于 地 形 复杂 的 地 区 ,每 
个 风向 呈现 的 粗糙 元 素 的 密度 和 分 布 不 同 ”””” 。 
对 于 均匀 ,平坦 的 下 垫 面 来 说 , 随 着 风向 的 改变 , 空 
气动 力学 粗糙 度 变化 幅度 较 平 组 ,这 是 由 于 影响 风 
浪 区 的 周 于 障碍 物 较 少 ,不 同 风向 上 风浪 区 内 粗糙 
元 差别 不 大 ,因此 均匀 下 垫 面 可 以 不 考虑 风向 的 影 
响 ; 而 对 于 非 均匀 下 垫 面 ,不 同 风向 造成 地 表 粗 糙 元 
对 滑 流 的 拖 搜 力 不 同 ,导致 不 同 风向 下 的 空气 动力 
学 粗糙 度 变 化 , 且 空 气动 力学 粗糙 度 随 风向 而 变化 
的 程度 也 不 同 ,在 特定 风向 范围 内 , 随 风 向 变化 幅度 
较 大 下。 因此 ,空气 动力 学 粗糙 度 随 风向 变化 是 
由 风浪 区 地 表 粗 糙 元 和 气流 共同 作用 形成 的 ,很 难 
通过 观测 获得 非 均 匀 下 执 面 的 空气 动力 学 粗 燃 度 。 
于 名 如 “基于 气象 站 数据 和 通 量 观 测 数 据 , 发 
现 同一 季节 的 空气 动力 学 粗糙 度 值 在 不 同 风 向 上 
表现 出 明显 的 差异 ;Yu 等 ”基于 黑河 流域 不 同 地 表 
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表 粗 糙 元 变化 更 敏感 ”。 

描述 庇 流 的 许多 微 气象 参数 与 植被 冠 层 结构 
相关 ,空气 动力 学 粗糙 度 与 平均 植被 高 度 存 在 函数 
关系 "0。 植 被 下 垫 面 的 空气 动力 学 粗糙 度 随 植 
被 生长 状况 出 现 季 节 性 变化 ,当地 表 植 被 高 度 逐 渐 
升 高 时 ,植株 对 于 气流 的 拖 遇 力 增 大 ,空气 动力 学 
粗糙 度 变 大 。Li 等 “利用 激光 雷达 和 成 像 光谱 
技术 估算 空气 学 粗 燃 度 , 发 现 空气 动力 学 粗 烽 度 的 
准确 性 很 大 程度 上 取决 于 佑 计 模 型 和 灌木 在 模型 
中 的 高 度 。 

植被 密度 和 盖 度 反映 了 植被 的 分 布 .成 密 程度 
以 及 植物 进行 光合 作用 面积 的 大 小 。 植 被 密度 对 
空气 动力 学 粗糙 度 有 很 强 的 影响 ,对 于 相同 高 度 但 
密度 不 同 的 植被 , 较 稀 玻 的 植被 的 空气 动力 学 粗糙 
度 较 低 ”。 研 究 认为 扩大 植被 盖 度 能 提高 地 表 空 
气动 力学 粗糙 度 “"” ,这 是 由 于 植被 区 的 流 场 存在 
绕 流 ,植株 密度 对 流 场 的 影响 较 大 , 比如 在 夏季 ,地 


格局 的 3 个 自动 气象 站 的 数据 ,分 析 了 不 同 季节 影 
响 空 气动 力学 粗糙 度 的 主要 气象 因素 ,发 现在 草原 
和 和 森林 山区 ,风向 是 主要 影响 因素 ,但 对 于 均匀 植 
被 下 垫 面 , 风 向 可 以 忽略 不 计 。 
3.1.4 大 气 稳定 性 空气 动力 学 粗糙 度 与 近 地 表 的 
大 气 层 结 有 关 , 大 气 层 结 一 般 用 稳定 度 来 度量 ,不 
同 的 大 气 稳定 度 影 响 满 流 的 产生 “>” 中 。 以 往 对 
于 空气 动力 学 粗糙 度 的 研究 是 基于 大 气 处 于 中 性 
稳定 的 条 件 下 ,但 在 通常 情况 下 ,这 种 理想 状态 难 
以 实现 。 不 同 大 气 稳定 条 件 下 的 空气 动力 学 粗粮 
度 差异 明显 ,于 名 召 " 发 现 相同 地 表 条 件 下 ,大 气 
处 于 稳定 条 件 下 的 空气 动力 学 粗糙 度 明 显 大 于 不 
稳定 条 件 , 且 不 稳定 大 气 条 件 下 空气 动力 学 粗糙 度 
的 离散 程度 较 小 。 大 气 稳定 性 与 昼夜 交替 紧密 相 
关 , 日 天 , 受 太阳 辐射 影响 ,由 于 植被 与 大 气 之 间 发 
生 消 流 交换 ,大 气 条 件 不 稳定 ,空气 动力 学 粗糙 度 
保持 低 水 平 ;夜间 ,大 气 注 流 作用 减弱 ,大气 处 于 较 
为 稳定 的 条 件 下 ,空气 动力 学 粗糙 度 显著 增加 ”。 
3.2 地 表 粗 糙 元 形态 特征 因素 

地 表 粗 糙 元 (如 植被 及 其 他 障碍 物 ) 的 结构 和 
空间 分 布 是 空气 动力 学 粗糙 度 的 影响 因素 ,植被 的 
形态 特征 包括 植被 高 度 .密度 和 盖 度 等 ,影响 地 表 
粗糙 状况 及 气流 拖 忠 力 大 小 ”™*“*。 研 究 指 出 ， 
空气 动力 学 粗糙 度 与 大 气 层 结 稳定 度 有 关 , 但 对 地 


表 植 被 盖 度 相对 较 大 ,植被 株数 较 多 ,所 以 粗糙 单 
元 分 布 密 度 相 对 较 大 ,空气 动力 学 粗糙 度 将 随 着 流 
场 阻力 的 增 大 而 增 大 。 

LAI 对 空气 动力 学 粗糙 度 的 影响 存在 规律 性 。 
在 植被 生长 初期 ,LAI 较 小 ,叶片 稀疏 有 旦 柔软 ,对 气 
流 的 阻力 较 小 ,空气 动力 学 粗糙 度 较 小 ; 随 着 植被 
不 断 生长 ,空气 动力 学 粗糙 度 随 LAI 增 大 而 增 大 , 达 
到 一 个 峰值 后 再 降低 。 原 因 是 当 植被 生长 达到 一 
定 浓密 程度 ,由 于 植被 叶片 过 于 密集 , 风 不 能 从 作 
物 内 部 穿 过 ,直接 从 植被 顶部 流 过 ,反而 降低 了 对 
气流 的 阻力 ,植被 的 表面 粗糙 度 降低 ,使 得 空气 动 
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遥感 技术 对 大 尺度 地 区 空气 动力 学 粗糙 度 参 
数 化 葛 定 了 基础 ,开辟 了 多 种 研究 方法 ,目前 遥感 
观测 系统 在 时 间 、 空 间 和 光谱 分 辨 率 上 正在 不 断 革 
新 升级 ,如 何 更 好 应 用 到 空气 动力 学 粗糙 度 领域 需 
要 综合 考虑 地 表 空 气动 力学 特性 。 由 于 多 角度 观 
测 可 以 观测 到 植被 对 太阳 光 的 不 均匀 散射 ,多 角度 
和 遥感 数据 在 描述 作物 冠 层 结构 方面 具有 不 可 比拟 
的 优势 。 因 此 , 随 着 遥感 技术 向 多 角度 观测 .高 时 
空 分 辨 率 发 展 ,遥感 反 演 的 空气 动力 学 粗糙 度数 据 
集 的 准确 性 得 以 提高 。 
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遥感 方法 在 空气 动力 学 粗糙 度 领 域 显 示 了 其 
特有 的 优势 ,但 使 用 遥感 空间 分 析 方 法 在 宏观 太 度 
上 估算 植被 下 垫 面 空气 动力 学 粗糙 度 时 ,参数 化 时 
会 产生 不 确定 性 ,原因 主要 在 于 未 充分 考虑 地 表 起 
TR .植被 变化 及 地 气 能 量 平衡 动态 。 人 研究 证 明 ,地 
表 起 伏 会 影响 地 表 粗 糙 元 分 布 和 地 表 通 量变 化 ,对 
地 表 粗 烙 度 产生 影响 ,因此 空气 动力 学 粗糙 度 的 研 
究 应 综合 考虑 气象 因素 .粗糙 元 形态 学 因素 以 及 地 
形 因素 ,分 析 空 气动 力学 粗糙 度 的 主要 驱动 因子 ， 
从 多 维 数据 中 提取 有 效 特征 信息 。 未 来 需要 充分 
发 挥 多 种 数据 源 合理 交叉 融合 的 优势 ,发 展 更 为 精 
细 的 长 时 序 空气 动力 学 粗糙 度 变化 的 数据 集 ,整合 
遥感 .站 点 观测 与 地 表 滑 流 模型 模式 ,深入 分 析 空 
气动 力学 粗糙 度 在 不 同时 空 太 度 的 变化 。 对 于 复 
杂 下 垫 面条 件 下 空气 动力 学 粗糙 度 的 计算 ,应 注意 
实测 数据 的 验证 和 标定 ,将 遥感 反 演 与 地 面 观 测 绪 
合 ,确保 模型 的 准确 性 。 目 前 ,对 空气 动力 学 粗粮 
度 的 佑 算 和 分 析 人 研究 主要 集中 在 植被 下 垫 面 ,未 来 
可 以 对 城市 等 其 他 下 垫 面 类 型 的 空气 动力 学 粗糙 
度 进行 进一步 拓展 人 研究。 
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Research progress on aerodynamic roughness 


LI Xinyu, WANG Jingpu, WANG Zhoulong 


(School of Resources and Environmental Engineering, Ludong University, Yantai 264025, Shandong, China) 


Abstract: Aerodynamic roughness is defined as the height at which the wind speed becomes zero under neutral 
and stable conditions. It is an important parameter for measuring the momentum and energy exchange between 
the underlying surface and atmosphere, and it is critical for investigating various surface processes and climate 
change. However, it has always been difficult to estimate aerodynamic roughness accurately at the regional scale, 
and there is no unified estimation model presently. Therefore, the parameterization of aerodynamic roughness is a 
topic worthy of further study. As a long-range monitoring method, remote sensing technology has the advantages 
of macroscopic and rapid acquisition of ground feature information, and its ability to achieve dynamic monitoring 
at the regional scale or a larger scale in estimating the aerodynamic roughness of vegetation-covered surfaces. 
Therefore, using remote sensing technology to estimate aerodynamic roughness has become a hot issue at home 
and abroad in recent years. In this study, the progress of research on aerodynamic roughness at home and abroad 
in recent years is systematically described. The estimation methods are divided into two categories: one is based 
on measured data, and the other is the remote sensing method, which is rapidly advancing. This study primarily 
introduces the method of estimating the aerodynamic roughness of the underlying surface of vegetation by remote 
sensing technology. Methods based on measured data include the canopy height fixed ratio method, field experi- 
ment method, and wind tunnel method; remote sensing methods include vegetation index, LIDAR, and multi- 
source remote sensing synergistic methods. In addition, the advantages and disadvantages of the different meth- 
ods are summarized at the end of each section. Finally, this study analyses the influence of meteorological factors 
and morphological characteristics of surface roughness elements on aerodynamic roughness and discusses the de- 
velopment trends and problems of remote sensing techniques in estimating aerodynamic roughness, aiming to pro- 
vide ideas for subsequent research on remote sensing monitoring of aerodynamic roughness. 


Key words: aerodynamic roughness; remote sensing; research status; influencing factors 


